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Abstract. In order to solve the problem of eternal acceleration, a model has been 
recently proposed including both a negative cosmological constant A and a scalar field 
evolving under the action of an exponential potential. We further explore this model by 
contrasting it against the Hubble diagram of Type la supernovae, the gas mass fraction 
in galaxy clusters and the acoustic peak and shift parameters. It turns out that the 
model is able to fit quite well this large dataset so that we conclude that a negative A 
is indeed allowed and could represent a viable mechanism to halt eternal acceleration. 
In order to avoid problems with theoretical motivations for both a negative A term and 
the scalar field, we reconstruct the gravity Lagrangian f(R) of a fourth order theory 
of gravity predicting the same dynamics (scale factor and Hubble parameter) as the 
starting model. We thus end up with a f(R) theory able to both fit the data and solve 
the problem of eternal acceleration without the need of unusual negative A and ad hoc 
scalar fields. 



PACS numbers: 04.50.+h, 98.80.-k, 98.80.Es 
1. Introduction 

Recent astrophysical data, from the Hubble diagram of Type la Supernovae (SNela) 
[U, [2] to the measurement of the anisotropy and polarization spectra of the cosmic 
microwave background radiation (CMBR) [3111] and large scale structure data [5], point 
out towards a picture of the universe unexpected only few years ago. According to this 
new scenario, the universe is spatially flat, undergoing a phase of accelerated expansion 
and dominated by a smoothly distributed negative pressure fluid dubbed dark energy. 
Although a lot of candidates has been proposed for this mysterious component [6], the 
concordance ACDM model [7] made out of cold dark matter (CDM) and the Einstein 
cosmological constant A stands as the best fit to a wide range of different astrophysical 
data [El [91 [in]. Notwithstanding its observational merits, the ACDM scenario is seriously 
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plagued by the well known coincidence and fine tuning problems which are the main 
motivations to look for alternative models. 

Although as a classical (general relativistic) field theory the concordance model is 
fairly simple, a universe presently dominated by a positive A term is on the contrary 
quite difficult to understand from a quantum theory viewpoint. Indeed, since matter and 
radiation energy density fades away as a~ 3 and a~ 4 respectively (with a the scale factor), 
the universe turns out to be asymptotically dominated by the vacuum energy accounted 
for by the positive A. As a consequence, the universe asymptotically enters a de 
Sitter phase with the scale factor growing exponentially and a never ending accelerated 
expansion. Let us then consider the definition of cosmic horizon as : 



with t and t en d the present day age of the universe and its age at the end of time (i.e., 
a{t e nd = 0). Note that, while for a closed universe t en d takes a finite value, this is not 
for open and spatially flat models (t en d — > oo). For a de Sitter like universe, S takes a 
finite value so that a cosmic horizon appears. It is worth stressing that this problem 
is not strictly related to A, but turns out to be a consequence of eternal acceleration. 
Should we replace A with a quintessence scalar field, the universe should still be eternally 
accelerated finally reaching a de Sitter phase and hence again a finite cosmic horizon. 

The presence of a cosmic horizon implies that it is not possible to define pure state 
vectors of quantum asymptotic states. As a consequence, in a de Sitter background 
spacetime, it is not possible to define a gauge invariant scattering matrix S. First 
studied in [TT], this problem is known as superexpansionary divergence in the context 
of quantum field theory (QFT) and is particularly troublesome for the formulation 
of strings theory. Indeed, by construction [12J, perturbative string theory is based 
on the well defined nature of scattering amplitudes of various excitations and hence 
fundamentally relies on the possibility to correctly define the S- matrix [13]. As such, 
formulating string theory in presence of a finite cosmic horizon is challenging. The so 
called Liouville string framework [TJ] represents an attempt to solve the problem in the 
case of non conformal backgrounds (including the de Sitter spacetime), but much work 
is still needed before a complete mathematical formulation is achieved. 

A radically different approach is, however, possible. Indeed, observations do tell us 
that the universe is spatially flat and accelerating, but this by no way implies that the 
acceleration should be eternal. This latter feature turns out only as a consequence of 
assuming a positive A terrrj§| to explain cosmic speed up. From an observational point 
of view, however, nothing prevents us to consider an effective dark energy fluid as a 
source of the accelerated expansion. We can then split its energy density as the sum of 
two components and then investigate the consequences on the future evolution of the 
universe. Such an approach has been recently investigated by some authors [T5l [16j [17] 
with interesting results. In particular, in [16] , some of us have presented a model whose 

§ A similar argument also holds for a quintessence scalar field whatever is the self interaction potential. 
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dark energy component is the sum of a negative cosmological constant and a quintessence 
scalar field evolving under the action of an exponential potential. As a result, although 
the model is presently accelerating, eternal acceleration disappears and the universe 
ends in a Big Crunch like singularity in a finite time. Motivated by these theoretical 
virtues, we further explore this model by the observational point of view in order to see 
whether a negative A is indeed compatible with the astrophysical data at hand. 

Notwithstanding the positive evidence of the model, it is also worth wondering 
whether an alternative theoretical derivation is possible. Actually, there are difficulties 
to reconcile a negative A term within the framework of Quantum Field Theory. In the 
context of the Standard Model of Elementary Particle Physics, Spontaneous Symmetry 
Breaking induces a large (negative) value of the cosmological constant: 

M 2 H 

md ~~8V2G~F () 
If we apply the current experimental bounds on the mass of the Higgs boson Mh and on 
the Fermi constant Gp, we obtain a (negative) value 55 orders of magnitude larger than 
the observed upper bound of the cosmological constant CC. There have been several 
attempts to solve this problem, without enough success [18J. A vacuum cosmological 
term A„ with (positive) opposite sign is often introduced, so that the physical observable 
cosmological constant results : 

A p h = A ind + Ay (3) 

The introduction of A„ is also dictated by the requirement of renormalizability of the 
massive theory, so the CC problem is rather the need of the extremely precise choice 
of the corresponding normalization condition A^ = in the very far infrared. It is 
very difficult to explain why the two terms on the right hand side of above equation 
should cancel each other with such accuracy at that point. On the other hand, as for 
all quintessence models, it is not totally clear where the scalar field comes from. Having 
such a large dominance by dark ingredients is, for sure, a disturbing situation so that it 
is worth to look for other potential explanations. 

Following the prescriptions of the Occam razor, it is desirable to have as less 
components of unknown origin as possible in a cosmological model. In order to overcome 
this problem, we resort to fourth order theories of gravity (also referred to in the following 
as f(R) theories). According to this approach, cosmic acceleration is the first signal of 
a breakdown of General Relativity so that the gravity Lagrangian has to be modified by 
replacing the scalar curvature R with a generic analytic function f(R) [HI [20j ED 122] . 
Moreover, it has been shown [23J that it is always possible to find out a f(R) theory 
reproducing the same dynamics (i.e., scale factor and Hubble parameter) as a given 
dark energy model. Using this method, we can reformulate the model above in terms of 
fourth order theories of gravity. The resulting f(R) will allow to both fit the data and 
solve the problem of eternal acceleration without resorting to a controversial negative 
A and an ad hoc scalar field. 
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The plan of the paper is as follows. We present the starting cosmological model 
in Sect. 2 where we generalize the approach in [16] including the radiation term. The 
numerical solution for the expansion rate for this model is the key ingredient for the 
likelihood analysis discussed in Sect. 3 where we contrast the model against a wide set of 
observational data thus being able to constrain its characteristic parameters. Although 
our model comprises two ingredients other than standard matter and radiation, it is 
worth stressing that it is dynamically equivalent to a model dominated by a single dark 
energy fluid whose equation of state is reconstructed in Sect. 4. The best fit model 
is, on the other hand, used as input for the procedure of reconstructing f(R) which is 
presented in Sect. 5 where we show that a negative A could result as part of an effective 
representation of a fourth order theory. Sect. 6 is finally devoted to conclusions. 

2. The model 

Looking at the impressive amount of papers addressing the problem of cosmic 
acceleration clearly shows that two leading candidates to the dark energy throne are 
the old cosmological consntat A and a scalar field <ft evolving under the influence of its 
self - interaction potential V((f>). While a A term adds a constant energy density p\ and 
a negative pressure = —p Lambda into the cosmic evolution equations, a scalar field 
partecipates to the dynamics through its energy density and pressure given as : 

' = \tf + v(4>) 

(4) 

p* = U 2 -v{4>) 

with a dot denoting the derivative with respect to the cosmic time t. The evolution of 
the scalar field is then governed by its Klein -Gordon equation: 

dV 

+ 3#0+ — = (5) 
d<p 

with H = a I a the Hubble parameter and a the scale factor (normalized to unity at the 
present time). It is easy to show that Eq.(jS]) is the same as the continuity equation for 
the scalar field energy density p^ under the hypothesis that (j) does not interact with the 
matter and the other ingredients of the cosmic pie. 

In the usual approach, one adds either a scalar field or a cosmological constant 
term to the field equations. However, since what we see is only the final effect of the 
dark energy components, in principle nothing prevents us to add more than one single 
component provided that the effective dark energy fluid coming out is able to explain 
the data at hand. Moreover, as we have hinted upon above, a single scalar field, while 
explaining cosmic speed up, leads to a problematic eternal acceleration. A possible way 
out of this problem has been proposed by some of us [16] through the introduction 
of a negative cosmological term. Motivated by those encouraging results, we therefore 
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consider a spatially flat universe filled by dust matter, radiation, scalar field and a 
(negative) cosmological constant term. The Friedmann equations thus read : 



H 2 



1 



PAI + Pr + PA + 7^<fi 2 + V{4>) 



2H + 3H 2 = - 



\pr -PA + \ft ~ V{4>) 



(6) 
(7) 



where we have used natural units with 8nG — c — 1. Using the continuity equation for 
Pm and p r and the definition of density parameter for the i-th fluid 



Pi{a 



Pi{a 



3H 2 



Eqs.flEJ) and (J7|) may be rewritten as: 



E 



ml \ alt 



». 2 



H n ^E = — Vt\ja ^ — 1Vt r a 4 — 



IHl \ alt 



(8) 



(9) 



with E = H/Hq, <j) = 4>/V3, V = V/3Hq, and we denote with a subscript the present 
day value of a quantity. It is convenient to introduce the dimensionless variable 

u = In (1 + z) = — In a 

with z = 1/a — 1 the redshift (having set a^ = 1). In terms of this variable, the 
Friedmann equations become : 



2 n M e 3u + n r e 4u + n A + v((j)) 



E 



lalE 2 



(1/2) (al^/alu)' 



^n M e 3u + 2n r e Au + Z -E 2 \ — 



,7\ 2 



2 rfw 

while the Klein- Gordon equation reads: 
,al 2 d) (lalE 2 x - 



E' 



alu- 



+ 



2 



3£r 



alu 



du 



dV 
+ -~ = . 



(10) 



(12) 



As well known, given an expression for the scalar field potential V(0), only two out of 
the three equations (jlOl . (fill) . (fl2|) are independent. It is convenient to insert Eqs.( fT0l) 
and ( TTTT) into Eq. (fl2l) to get a single equation governing the evolution of the scalar field : 



,d 2 ({ 
du- 



+ 



Q M e 3u + Q r e 4u + 3Q A + 3V(<f>) 



du 







(13) 
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with E 2 (u) given by Eq. fllOp . In order to solve (numerically) this equation thus 
determining <p{u) and then E(u) through Eq. ffTOj) . one has to set two initial conditions. 
The first one can be trivially obtained by evaluating Eq. ( TTOT) in u = to get : 



u=0 



with Vq = V((f> ). Given the shape of the potential V(<p), the above relation can be 
inverted to get the present day value of the scalar field. A second initial condition can 
be obtained by first considering the equation of state (hereafter EoS) of the field given 
as : 

n E*(dHdu) 2 -2V 

™* = — = ^ • ( 14 ) 

P<t> E 2 [d(f)/du) + 2V 

Evaluating this at the present day and using the above relation for </> , we finally get 
the initial conditions : 

V = ^(1 + ioo) , (15) 

g))' =*(!-«.) (16) 

with wo = w^iu = 0). Note that, since the left hand side of Eqs.( |T5|) and ( |T6l) are 
positive definite, we get the constraint — 1 < wo < 1 which is always verified for any 
ordinary scalar field whatever the potential V(<j>) is. 

2.1. The deceleration parameter 

In order to separate models in accelerating and decelerating ones, one has to compute 
the present day value of the deceleration parameter defined as : 

da H 1 dE 

q{u) = -^ = - 1 -lP = - 1 + Elu- 



Using Eqs. fllOp and fllip . we can therefore write the deceleration parameter as: 



~\ 2" 

17) 



x 



\du J 

(3/2)tt M e 3u + 2VL r e 4u + (3/2) (d&fdu 



Evaluating this relation in u = and using Eq.(Tl6l). the present day value then reads: 



1 3 
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so that, in order to have accelerating models (go < 0), we have to set the constraint : 

wa £ gA - (i + Ml . (19) 

We can also solve Eq. (ITH1) with respect to u; to get : 

(l-2g ) + ^ r -3^ A 

W0 = 3^ (20) 

so that the initial conditions now rewrite : 

V Q = (3fi + Vl r - 3fi A + 1 - 2q )/Q , (21) 



du 



\ =3tt0-tt r + 3tt A - l + 2q . (22) 

u=0 

Since the lhs of these relations are positive quantities, we get the following constraints 
on the present day deceleration parameter : 

-(i + n r ) + + n A ) < qo < -(i + o r ) + -(^ - q a ) . (23) 

For accelerating models, the lower limit in Eq. f[2"3"j) must be negative so that we get : 

^A<~(l + ^r)-£V (24) 

Note that for models with a negative cosmological constant this relation is always 
satisfied thus meaning that for such a choice it is possible to work out accelerating 
solutions of the field equations. 

2.2. The exponential potential 

The above discussion holds whatever is the scalar field potential V(<f)), but in order to 
actually solve the cosmic equations we have to definitively assign an analytical expression 
for this quantity. Following [241 [25] , we consider an exponential potential : 

V((f>) = B 2 exp(-a<f)) (25) 

with B 2 a generic constant and cr 2 = 12ttG/c 2 = 3/2. As shown in [26], such a choice 
leads to a cosmological model (made out of the above scalar field and dust matter) in 
good agreement with a large set of astrophysical data. 

Although being interesting on both theoretical and observational grounds, such a 
model is however affected by the problem of eternal acceleration. As can be easily 
understood, since the matter energy density decreases as a~ 3 (with a the scale factor 
normalized to be unity at the present day), the universe becomes soon scalar field 
dominated and the accelerating expansion never ends. In order to avoid this problem, 
Cardenas et al. [IS] have added a third ingredient to the cosmic pie, namely a negative 
cosmological constant. The Nother symmetry approach [27] makes it possible to find a 
convenient change of variables in such a way that exact solutions are found. 
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Although dealing with an analytical expression is welcome, the result in [TO] only 
holds if the radiation term is neglected. While this is not a problem when considering 
the late universe, introducing radiation drastically changes the structure of the cosmic 
equations so that the Nother symmetry approach does not apply anymore. As a 
consequence, the analytical solution may no more be used and a numerical analysis 
is needed. To this end, we insert Eq.fl25l) into Eq. (fT5]) and solve with respect to (fi : 




l^ka-^o) 



V 9 ln ~ W °^ 



(26) 



where, in the second equality, we have arbitrarily set B 2 = (3/2)Hq. Note that such a 
choice has no effect on the dynamics since it is a simple rescaling of the scalar field <fi 
which does not influence any physically interesting quantity. With such a choice, the 
potential and its derivative then read : 



v(4>) 



dV 




(27) 



Summarising, in order to determine the cosmic dynamics (i.e., the scale factor and the 
Hubble parameter) for such a model with both a negative cosmological constant A and 
a scalar field (ft with an exponential potential added to the standard dust matter and 
radiation terms, we have to numerically integrate Eq. (fl3|) with the potential given by 
Eq.(l2"TI) and the initial conditions (I2"6"j) and (TTBT) . The solution for 0(w) thus obtained 
may then be inserted into the potential and then into Eq. (1101) to determine the Hubble 
parameter which can be further integrated to get the scale factor. In order to perform 
such a scheme, one has to set the values of up to four parameters, namely the dust and 
radiation density parameters Qm and fl r , and the present day values (Q^jWo) of the 
scalar field energu density and EoS respectively. In particular, one could use Eq.(fT9l) to 
choose Wq in such a way that the resulting model will be today accelerating. 



3. Matching with the data 

Notwithstanding how well motivated it is, a whatever model must be able to reproduce 
what is observed. This is particularly true for the model we are considering because 
the presence of a negative cosmological constant introduces a positive pressure term 
potentially inhibiting the cosmic speed up. Moreover, contrasting the model against the 
data offers also the possibility to constrain its characteristic parameters and estimate 
other derived interesting quantities, such as go, the transition redshift zt and the age 
of the universe to. Motivated by these considerations, we will therefore fit our model 
to the dataset described below parametrizing the model itself with the matter density 
Qm, the scalar field quantities (Q^,wq) and the dimensionless Hubble constant h (i.e., 
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H i n units of 100 km/s/Mpc), while we will set the radiation density parameter as 
Q r = 10~ 4 ' 3 as in [28] from a median of different values reported in literature. 

3.1. The method and the data 



In order to constrain the model parameters, we maximize the following likelihood : 



£ oc exp 



x 2 (p) 



(2f 



where p = (Qjvf; ^4>i w o, h) denotes the set of model parameters and the pseudo-x 2 
merit function reads : 



N 



x 2 (p) = E 



i=l 



fi th (z t ,p) - ^ bs 



+ E 



N r fth I ~ n \ _ -fobs 
J gasK^ii trj J gas,i 



i=l L 



-i 2 



+ 



A{p) - 0.474 
0.017 



+ 



K(p) - 1.70 
O03 



'h- 0.72 s 
0.08 



(29) 



Let us discuss briefly the different terms entering Eq. ff29|) . In the first one, we consider 
the distance modulus fj, — m — M, i.e. the difference between the apparent and absolute 
magnitude of an object at redshift z. This is given as : 



m 



M = 25 + 5 log D L (z) 



with D L (z) the luminosity distance (in Mpc) defined as: 

dz' 



D L (z) = ^(l + z) r 

-fin JO 



E z' 



(30) 



(31) 



As input data, we use the SNela sample assembled in [29j by putting on a common scale 
the data recently released from the SNLS [30] and ESSENCE collaborations [31J and the 
higher redshift SNela observed with HST in the GOODS survey [32]. As well known, 
the SNela Hubble diagram is unable to determine the Hubble constant Ho since this 
quantity is degenerate with the (unknown) absolute magnitude M. As such, Davis et 
al. have put all the SNela on the same distance scale using h = 0.656. In our analysis, 
however, we will leave h as a free parameter so that Eq. fl30l must be rewritten as : 

/i(z) =25 + 5 log + 5 log f ^7— 1 + 5 log d L (z) (32) 



Tjf id 



h 



with Hl %d 



100/i^ id = 65.6 km/s/Mpc and d^ = Dl/(c/Hq) the Hubble free luminosity 
distance. Note that the value of h determined by this fit will not be fully reliable because 
of the degeneracy hinted above. However, such a problem will not affect the estimates 
of the other parameters since we will marginalize over h in the analysis of the results. 

The second term in Eq.( f29i) relies on the gas mass fraction in galaxy clusters. We 
briefly outline here the method referring the interested reader to the literature for further 
details [33J, EH [35] . Both theoretical arguments and numerical simulations predict that 
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the baryonic mass fraction in the largest relaxed galaxy clusters should be invariant with 
the redshift (see, e.g., Ref. [36]). However, this will only appear to be the case when 
the reference cosmology in making the baryonic mass fraction measurements matches 
the true underlying cosmology. From the observational point of view, it is worth noting 
that the baryonic content in galaxy clusters is dominated by the hot X-ray emitting 
intra - cluster gas so that what is actually measured is the gas mass fraction f gas and it is 
this quantity that should be invariant with the redshift within the caveat quoted above. 
Moreover, it is expected that the baryonic fraction in clusters equals the universal ratio 
Q b /Q M so that f gas should indeed be given by bxfl b /Q M where the multiplicative factor 
b is motivated by simulations that suggest that the gas fraction is slightly lower than 
the universal ratio because of processes that convert part of the gas into stars or eject 
it outside the cluster. 

Following Ref. [37], we adopt the concordance ACDM model (i.e., with Qm = 0.3, 
= 0.7, h = 0.7) as reference cosmology in making the measurements so that the 
theoretical expectation for the apparent variation of f gas with the redshift is [37] : 

3/2 

(33) 



K 7 A{z)b(z) n b 

Jgas\Z) - 



i + s(z) n 



M 



D\ CDM {z) 



Da(z,p) 



Some words are needed to explain the meaning of the different terms entering the 
above equations. First, the two functions b(z) and s(z) take into account variations 
of the ratio Q b /flM due to gastrophysics and star formation respectively. In a first 
approximation,, they are described as linear functions of the redshift. However, in 
order to not increase the number of parameters, we will set them to constant values 
taking b(z) = b = 0.83 and s(z) = s = 0.16(/i/0.70) 1/2 in agreement with [37J. While 
K and 7 are normalizing factors correcting for measurement related problems and can 
be put to the constant values 1.0 and 1.05 respectively, A(z) actually plays a more 
important role. To understantd its origin, it is worth remembering that f gas is typically 
measured at a given fraction of the cluster virial radius. However, on the sky distances 
are measured in angular rather than physical units with the conversion depending on the 
assumed cosmological model. Since the reference model is different from the actual one, 
a correction term must be included to account for the change in the radius. Considering 
the scaling of the different quantities involved in the measurement process, it is possible 
to show that the correction term may be approximated as 



A{z) 



E(z)D A (z,p) 



E ACDM (z)D JVCDM (z) 



(34) 



with E ACDM (z) = [Fl M (l + zf + Q A } 1/2 and 77 = 0.214. Note that, actually, because of 
the small 77 value, this term does not play a significant role in the fitting process, while 
a key ingredient is the last one in Eq. fl3"3"|) given by the ratio of the angular diameter 
distances D^(z) = D L (z)/(l + z) 2 between the ACDM and the model to be tested. It is 
worth noting that f gaa (z) depends not only on the integrated Hubble parameter, but also 
explicitly on the baryon and total matter density parameters Q b and Q&f. In particular, 
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baryogenesis calculations contrasted to the observed abundances of primordial elements 
puts a severe constraint on the physical baryon density Uj, = Vt^h 2 . Using this method, 
Kirkman et al. [38 1 have determined: 



u b = 0.0214±0.0020 . 

A rigorous analysis should be done letting uj, as a free parameter eventually including the 
above estimates as a prior. However, in order to not increase the number of parameters 
to be determined, we will set = 0.0214 neglecting the small error. 

The third term in the definition of x 2 takes into account the measurement of 
the baryonic acoustic oscillation (BAO) peak in the large scale correlation function 
at 100 h~ x Mpc separation detected by Eisenstein et al. [39] using a sample of 46748 
luminous red galaxies (LRG) selected from the SDSS Main Sample [ID]. Actually, rather 
than the position of acoustic peak itself, a closely related quantity is better constrained 
from these data, namely the acoustic peak parameter defined as [39J : 

n i/3 

(35) 



_ v"a7 
zlrg 



Z LRG 2 f 

-V {zlrg) 



E{z LRG y 

with zlrg — 0.35 the effective redshift of the LRG sample, and we have introduced the 
dimensionless coordinate distance y(z) defined as: 

Viz) = f ~~ppT~T~ \ • (36) 
Jo E{z,p) 

As it is clear, the A parameter depends not only on the dimensionless coordinate distance 
(and thus on the integrated expansion rate), but also on Qm and E(z) explicitly which 
removes some of the degeneracies intrinsic in distance fitting methods. Therefore, it is 
particularly interesting to include A as a further constraint on the model parameters 
using its measured value [3S] : 

A=OAm (o£) ±0Mr 

with n s the spectral index of the primordial density perturbations. For n s = 0.95 
as determined by the WMAP 3rd year analysis |41| . we get A = 0.474±0.017 as set in 
Eq. (|29l) . A caveat is in order here. As discussed in [39] , the measurement of the position 
of the BAO peak from the correlation function relies somewhat on having assumed the 
ACDM model to convert angular distances in physical distances and in the computation 
of the reference (smoothed) power spectrum. As a consequence, one should use with 
caution the above value for A when using a different cosmological model and rather 
directly fitting the measured correlation function. However, it is expected that the 
position of the BAO peak does not change too much in a different model since it is 
measured at a relatively low redshift. Since a full analysis will require the solution 
of the perturbation equations for our model (which is outside our aims here), we will 
follow the common practice in literature neglecting this problem and directly using the 
A parameter as an observational constraint included in the \ 2 merit function. 
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Let us consider the fourth term in Eq. (1291) which relies on the shift parameter [12] : 



with zls the redshift of the last scattering surface which we compute using the 
approximation given in [13]. Using the WMAP3 data, Wang & Mukherjee jB] have 
determined 1Z = 1.70±0.03 in very good agreement with what is expected for the 
concordance ACDM model. Note that such a result may argue in favour of models with 
a negligible dark energy component at high redshift, but we defer this discussion to the 
later analysis of the results. 

Finally, the term depending on h in Eq.( l29l) is only a Gaussian prior on this 
quantity obtained by considering the model independent estimate of the Hubble constant 
recovered by the HST Key project [45] . This collaboration have measured Hq using 
a wide set of different local distance calibrators thus ending up with a value which 
is claimed to fully take into account any possible systematic error inherent to the 
peculiarities of each single method. Comfortably, such a measurement turns out to 
be in good agreement with other (less precise) methods relying on different physics and 
distance scales as the time delay in lensed quasars [16] and the Sunyaev- Zel'dovich 
effect (17] in galaxy clusters. 

In order to maximize the likelihood function £(p), we should compute it over a very 
fine grid in the four dimensional parameter space (0^, O^,, w , h) and then interpolate 
the results for values falling in between two grid points. This is quite time consuming 
for a very fine grid so that we resort to a Monte Carlo Markov Chain (MCMC) method 
running three chains with 30000 points each and assessing the convergence using the 
Gelman-Rubin test (with \R — 1| = 0.1). It is worth stressing that, in order to be sure 
that the MCMC runs into the region with negative A, it is better to reparametrize the 
model using Q\ instead of fl^ as model parameter using the obvious relation : 



Therefore, in the following, we will use (flM, ^a, w , h) as the parameter space to be 
explored by the MCMC code. After cutting the burn -in period, the final coadded chain 
contains ~ 88000 points thus guaranteeing an efficient coverage of the interesting region 
of the parameter space allowing us to compute the marginalized likelihood functions for 
each parameter pi : 



which is then normalized at unity at maximum. Under the Bayesian framework, the best 
estimated for the parameter pi is given by the median of the marginalized likelihood, 
while the 68 and 95% confidence ranges are given as (xi ai yi a ) and (x 2(J ,?/2<t) with x ia 
and z/ia computed by solving respectively the equations : 




(37) 




1 — f^M — — • 




(38) 
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Table 1. Best fit (bf) and median (med) values and la and 2a ranges of the parameters 
(0m> ^Aj wq, h, 0,$) as obtained from the likelihood analysis. 



Par 


bf 


med 


la 


2a 


VLm 


0.283 


0.307 


(0.272,0.352) 


(0.246,0.410) 




-0.072 


-0.298 


(-0.54,-0.11) 


(-0.92,-0.02) 


w 


-0.72 


-0.67 


(-0.74,-0.60) 


(-0.79,-0.53) 


h 


0.632 


0.620 


(0.588,0.654) 


(0.554,0.692) 




0.789 


0.989 


(0.799,1.226) 


(0.700, 1.574) 



/Pi, max fPi,max 
C Pi (pi)dpi = Si / C Pi (pi)dpi , 

with 5i = (1 — 0.68)/2 for i = 1 and (1 — 0.95)/2 for i = 2 and (pi, m i n ,Pi,max) the lower 
and upper bounds chosen for the parameter p^. For the model we are considering, we 
conservatively let VLm range between 0.15 and 0.45, while the range for h is (0.45, 0.85). 
While the upper bound for Q\ is dictated by our constraint A < 0, choosing an upper 
limit is a more complicated issue. Indeed, the scalar field and the negative cosmological 
constant may be incorporated in a single dark energy fluid (as we will see later) so 
that there is a degeneracy in balancing the two individual components. We therefore 
arbitrarily set Qa > —1-5 thus allowing values of quite larger than 1. Finally, we 
cut the physically acceptable range (—1,1) for u> to become (—1,0) which is still a 
conservative choice. Indeed, we could also use Eq. (|T9|) to set an upper bound for w Q 
ensuring go < 0. However, we prefer to be fully open minded not forcing the chains to 
explore a priori acccelerating models only. 

3.2. Results 

Best fit model parameters, median values and 1 and 2a ranges for the parameters 
(fiM> ^a> w , h, are reported in Table 1, while Figs.[T]and[2]shows how well our best 
fit model reproduce the data on the SNela Hubble diagram and gas mass fraction. 

Let us first discuss some general points. First, we note that the best fit model is 
in quite good agreement with both the SNela and gas data. Indeed, the x 2 values are 
respectively 206 and 48 to be contrasted witth the number of datapoints, being 192 and 
42 respectively. Less good, but still in satisfactory agreement with the observed ones, 
are the values for the acoustic peak and shift parameters being : 

A = 0.45 , K = 1.67 . 

Motivated by these results, we can therefore safely conclude that including a negative A 
leads to a model still in agreement with the data so that this approach to halting eternal 
acceleration is a viable one from an observational point of view. In a sense, this is not 
surprising given that the net effect of the scalar field and negative A is to provide a dark 
energy fluid with negative pressure dominating the energy budget. Neveretheless, it is 
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Figure 1. Best fit curve superimposed to the data on the SNela Hubble diagram. 
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Figure 2. Best fit curve superimposed to the data on the gas mass fraction. Note that 
the theoretical curve plots indeed fg as (z)x(h/0.7) 1 - 5 with h set to its best fit value. 



worth stressing that, while the one fluid model has to resort to a peculiar EoS, the two 
fluids scenario relies on quite simple ingredients such as the exponential potential and 
well motivated cosmo logical constant. 

Comparing the best fit and median values of each parameter, it is soon clear that 
they are well different. This is, however, not a peculiarity of the model or the analysis, 
but rather a typical consequence of the non Gaussian shape of the likelihood. Since the 
full likelihood £(p is not the product of the marginalized likelihoods C Pi (jpi), the best 
fit values pbf are not given by the set of parameters that maximize each single C Pi {pi). 
Moreover, the marginalized likelihoods are markedly asymmetric so that the median 
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values may be quite different from the maximum values. This is particularly the case 
for Q\ because of the constraint imposed by hand on the negativeness of A. It is worth 
stressing, however, that all these caveats are common to the Bayesian approach to every 
data fitting problem so that they have not to be considered as model shortcomings. 

A successful fit to the data could, however, be meaningless if the estimated 
parameters take unrealistic values. It is therefore worth discussing how our constraints 
compare with other previous estimates. Let us first consider, therefore, the matter 
density parameter Qm- It is remarkable that the values in Table I are larger than the 
WMAP3 best fit model [H], being Q M = 0.268±0.018. Taken at face values, the two 
results are markedly different, even if they well overlap at la. It is worth stressing, 
however, that the WMAP3 result reported above also depends on the cosmological 
model fitted to the data. For instance, adding the running of the spectral index to the 
fiducial ACDM scenario shifts upward the estimated Qm now giving Qm = 0.282±0.020 
in better agreement with our constraint. Moreover, our estimate also overlaps well 
with Qm = 0.27±0.04 from fitting the Chevallier- Polarski- Linder (CPL) model 
w(z) = w + W\zj{\ + z) [48J to a dataset similar to our own [29J, but not inlcuding 
the gas data. As such, we are confident that the higher than usual matter content of 
our model is perfectly reasonable and do not consider this as a worrisome shortcoming. 
It is worth noting that a higher value for Qm is somewhat expected. Indeed, in order 
to compensate for the attractive effect of the negative cosmological constant, a higher 
than usual contribute from the scalar field is needed. But such a large Q^ should drive 
the universe towards a kind of superacceleration that needs to be compensated for by 
a larger matter content. The balance among these contrasting effects then induces 
a pushes Qm towards larger values. As a final remark, we also note that the 1 and 
2a ranges are quite larger, while typically Qm is better constrained in literature. On 
one hand, this is a consequence of the degeneracy hinted above among three density 
parameters (Qm,Qa,Q<i>) rather than the only two, (Qm,Qde), typically entering dark 
energy models. On the other hand, we are not using here either the CMBR spectrum 
which probes the very high redshift (where Qm is the dominant term) or the LSS data 
(which are sensitive to the matter content). It is therefore not surprising that the error 
we achieve on Qm is so large. 

While it is meaningless to compare the value of to the other estimates in 
literature given that we are considering a negative rather than positive A, one could 
naively think that wq should be compared to other determinations. However, it is easy 
to understand that this is not the case. Referring to the CPL parametrization as a 
prototype for phenomenological EoS and using a similar dataset as our own, Davis et 
al. [29J have obtained Wq = — 1.1 jlg g. Compared with the values in Table 1, we find 
only a marginal agreement at 2a level. However, the CPL parametrization assumes 
a single dark energy fluid, while our model comprises two different fluids so that the 
two results are not comparable. Qualitatively, the smaller is w (in absolute value), 
the smaller is the driving force speeding up the universe expansion so that one could 
be surprised that our model fits so well the data. Actually, what really matters is not 
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the EoS, but rather than the pressure p. Considering its present day value, we have 
P(f,(z = 0) = WoQ^pcru = —0.57p crit which compares qualitatively well to the CPL value 
Pcpl = (1 — ^M)woPcrit = —0.37p cr i t . Moreover, one should also take into account the 
different h value which leads to different values of p cr it- Considering that the difference 
somewhat fades away when integrating the Hubble parameter to get the luminosity 
distance, it is therefore not surprising that the Hubble diagram is almost the same 
betweem the one fluid phenomenological CPL model and our two fluids scenario. 

We do not discuss here the constraints on h reported in Table 1 since, as we have 
yet said above, this quantity is degenerate with the SNela absolute magnitude. Note 
that h does not enter either in the acoustic peak nor the shift parameters so that they 
do not help in breaking this degeneracy. On the other hand, the gas mass fraction 
fgas(z) weakly depends on h since it enters the function A(z) defined in ( 1341) trhough 
Da(z). It is easy to check that A(z) oc (h/h\cDM)~ v so that for r\ = 0.214, we get 
indeed a negligible effect. With all these caveats in mind, we do not care much about 
our findings on h which turns out to be lower than the h ~ 0.72 preferred by the 
local measuremnts [15] and CMBR data [H]. Motivated by these considerations, we 
have therefore marginalized over h in the discussion of the results, although it is worth 
stressing that the difference between h and h\cDM = 0.7 has to be taken into account 
when comparing h dependent quantities between the two models. 

As a byproduct of the MCMC code, we can also estimate some interesting derived 
quantities to be contrasted again measurements in order to get some more tests on 
the viability of the model. To this aim, given a quantity g(p), we compute its value 
for the points belonging to the MCMC and estimate its median value and 1 and 2a 
ranges as described above. In order to both speed up the computation and reduce the 
correlation between neighbouring points, we have thinned the chains extracting one over 
eight points thus ending up with a sample of ~ 11000 values which is enough to get 
reliable estimates. 

First, we discuss the results for the acoustic peak parameter whose median value 
and 1 and 2a ranges read : 

.4 = 0.460 , A E (0.434, 0.490) , .4 G (0.407, 0.523) . 

Pending the question about the validity of its measured value in models other than the 
ACDM one, we can therefore conclude that A predicted by our scenario is in excellent 
agreement with what is inferred from the correlation function of LRG galaxies. It is 
nevertheless worth extending this analysis by computing the power spectrum and hence 
the correlation function predicted by our model, although this is outside our aims here. 
As a second interesting quantity, we consider the shift parameter for which we find : 

ft = 1.38 , 1Z E (1.21, 1.53) , 1Z E (1.06, 1.65) , 

these values being respectively the median and 1 and 2a ranges. A comparison with 
the measured value, 71 = 1.70±0.03, should make us conclude that the model is not 
able to fit this quantity. While this is somewhat expected on the basis of the equivalent 
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dark energy model described in the next section, a more subtle issue has to be taken 
into account. Looking at Eq. fl371) and considering the constraints on Qm in Table I, 
it is easy to conclude that the problem with 1Z originates from i/(zls) < Vacdm(zls)- 
We have checked that this is indeed the case comparing y(z) for our best fit model 
with the ACDM curve. Actually, y(z) deviates from yACDM {%) more and more as z 
increases. Although worrisome, the smaller value of TZ should not be considered as a 
strong motivation against our model. As explained in JUJJ [50], TZ is an approximation 
for the acoustic scale £ a which is defined as : 



with c s the sound speed. The shift parameter is obtained by £ a using a different 

— 1/2 

normalization and approximating the denominator as Vt M . This is correct only if 
the dark energy fades away with z becoming subdominant at the last scattering. While 
this is true for the ACDM and standard quintessence models, this is not for our model 
as we will see in the next section. One should therefore compute £ a for our case and 
compare with the estimated i a = 303.6^{;2 [5Dj - Unfortunately, we are unable to perform 
such a computation here because the numerical solution for <fi(z) and hence a(t) becomes 
seriously unstable for large z. As a consequence, the integral in Eq. (l3~9l) turns out to be 
unreliable so that we can not discuss anymore this point leaving it for further work. 

As a final issue, we also consider the age of the universe. Proceeding as before, we 
get (in Gyr) : 



with the same meaning as before for the values reported. Compared with t = 
13.75±0.15 Gyr estimated for the concordance ACDM model from a combined analysis 
of the three -year WMAP, SNela and BAO data [41] . our result is in very good 
agreement. It is nevertheless worth stressing that the WMAP3 constraint is strongly 
model dependent so that one should better compare with the model independent 
measurements. However, Krauss & Chaboyer [51] have estimated to = 12. 61^6 Gyr 
from globular clusters, while a similar value, to = 12.5±3.5 Gyr, has been determined 
from nucleochronology [52]. Given the larger error bars, it is not surprising that our 
estimate is in very good agreement with these values too thus giving a further strong 
evidence in favour of our scenario. 

4. Effective dark energy 

The model we are considering is made out of two non standard fluids, namely a negative 
cosmological constant and a scalar field. From a phenomenological point of view, 
however, we can work out a model comprising a single dark energy fluid giving rise 
to the same expansion rate as the present one. Ies EoS may be evaluated as : 



a 



?r(l + z ls )D a (zls) 
Iq 1+Zls c s da/(aa) 



(39) 



to = 13.9 , t e (13.0,14.9) , t G (12.2,15.9) , 



l+w eff {u) 



(2/3)d\nE(u)/du- Q M e 3u + |fi r e 4 " E~ 2 {u) 



(40) 



1 - (n M e 3u + n r e 4u )E- 2 (u) 
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Figure 3. Scalar field (solid line) and effective dark energy (dashed line) EoS as 
function of z for the best fit model. 

so that the energy density may be then computed as : 

Pde 0) I Peru = ^Di?exp ^3 — (41) 

Fig. [3] shows the reconstructed effective dark energy EoS for our model setting the 
parameters (Qm,^a,wq) to the best fit value determined by the MCMC code. There 
are some interesting issues than can be drawn from this plot. First, we note that w e ff(z) 
and w ( j ) (z) track each other being essentially equa{jj] both in the late (u < 2) and early 
(u > 4) universe. This is not surprising since, for the best fit values we are using, the 
negative cosmological constant is quite small. However, a similar behaviour is obtained 
also for other values of the model parameters so that this is a rather common feature. As 
such, we will not investigate in detail how the shape of w e ff(z) depends on (Qm, ^<^>, wq), 
but our discussion below for the best fit model is nevertheless quite general. 

As it is apparent, w e ff(z) cannot be fitted by the most common phenomeno logical 
EoS parametrization such as the CPL one. Indeed, while for the CPL ansatz dw/dz = 
(1 + z)~ 2 has a monotonic shape, in our model dw e ff/dz is not at all monotonic and 
actually also changes its sign. For the best fit model shown in the plot, w e ff(z) first 
decreases down to the almost ACDM value w e jf — —0.98 for z ~ 4.0 and then starts 
increasing crossing the dust value w e ff = at z ~ 15 to finally stays constant to the 
stiff matter value w e ff = 1 for z > 210. Varying (Q M ,Q^,w ) changes these values, 
but not the shape of u> e // which always asymptotes to 1. Notwithstanding this unusual 
feature, the present day values of w e ff and dw e ff/dz are well within the 2a contours 

|| Actually, in Fig. 02 the two curves for w ( j,(z) and w e ff(z) are almost perfectly superimposed so that 
they can be hardly distinguished. However, we have checked that the larger in (in abslute value), 
the larger is the difference between the two EoS. 
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one can obtain when fitting the CPL formula to the same dataset we have used here. 
This is not unexpected given that, for z < 1, Dl(z) for a linear EoS closely matches the 
luminosity distance- redshift relation predicted by our model. 

As an alternative approach, one could resort to a model independent reconstruction 
of the EoS (see, e.g., [53] and references therein). Comparing u> e // in Fig. [3] with Fig. 
3 (right panel) in [53] shows us that the effective EoS significantly departs from what 
can be reconstructed by the data. It is nevertheless worth noting that in both cases 
the EoS asymptotes w = 1 even if the transition is still more rapid in the case of 
the reconstructed EoS because of the lack of the initially decreasing phase we have 
in our model. Although troublesome, we are confident that the discrepancy between 
the reconstructed EoS and our w e ff(z) should not be considered as a strong evidence 
against the present model. Indeed, data do not tell anything on w(z) directly, but only 
on the luminosity (or angular diameter) distance as function of the redshift. As such, 
one should not be worry about not reproducing a given w(z), but the aim of the model 
must only be to fit the data with physically meaningful values of the model parameters. 
This task is successfully accomplished by our model so that we still retain it as a valid 
alternative to a single dark energy fluid scenario. 

A potentially more worrisome problem is represented by the behaviour of the density 
parameter £Ide(z) when compared to the matter one VLm{z). As it is apparent from 
Fig.Hl Q DE (z) closely follows f^(z) being only slightly larger because of the need for 
compensating the negative cosmo logical consntat term. From this plot, it is clear that 
the usual sequence dark energy - matter - radiation dominions is not realized in this model 
since the matter era is followed by another dark energy dominated epoch. As such, 
the model could be troublesome when dealing with the nucleosynthesis which takes 
place during the radiation era. Indeed, it is usually claimed that, in order the BBN to 
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be efficient, the constraints Mev) < 0.1 has to be imposed which is strikingly 

violated by our effective dark energy model. It is, however, worth noting that such a 
constraint usually refers to a negative pressure fluid, while w e ff is positive for our model. 
Exploring whether BBN can take place in such a background is outside our aim here, 
but we will return later on this topic in the conclusions. 

While the BBN epoch is far away in a redshift range where our model could not 
be trusted anymore (for reasons explained later), structure formation typically probes 
a nearer epoch of the universe evolution. In the usual scenario, the growth of structure 
takes place during the matter dominated era well after the universe has left the radiation 
epoch. However, in our case, the radiation dominated era has been replaced by a dark 
energy dominated one and, moreover, the matter term is the leading one for a shorter 
than usual period. For instance, in the concordance ACDM scenario, matter drives 
universe expansion from z ~ z eq ~ 3000, while, for the best fit model in Fig.HJ the 
matter epoch starts later at z — 157. Studying in detail the growth of structure and its 
impact on observable quantities such as the galaxy clusters mass function and abundance 
is outside our aims here and will be deeply investigated in a forthcoming publication. 
However, we can anticipate that this is a subtle issue. On the one hand, denoting with 
S the matter density fluctuations, one usually has : 

5 + 2H5 - AttGpm = (42) 

so that the dark energy only enters through its effect on the Hubble parameter H(z). As 
we have checked, AH/H = 1 — H/H^ C dm decreases from ~ 10% to ~ for z ranging 
between the beginning [z ~ 150) and the end [z ~ 0.45) of the matter dominated era, 
taking (Qm,^a, h) = (0.3,0.7,0.7) for the concordance ACDM model. As such, we do 
expect that the growth index / = din 5 /din a is almost the same as the ACDM one so 
guaranteeing a good agreement with this kind of data. However, Eq.( )42l) is recovered 
assuming that dark energy does not cluster. While this is reasonable for a negative 
pressure fluid, this could not be the case for our model since its EoS increases from the 
negative value w e ff(z = 0.45) = —0.84 to w e ff(z = 150) ~ 0.99 during the matter era. 
As such, the dark energy sound speed c 2 s = dpnE/dpcE crosses the dust value c s = 0. 
As such, it is likely that dark energy can partecipate to the collapse leading to structure 
formation thus invalidating Eq. (l42p and claiming for a coupled set of equations. 

As a final issue, we would like to comment on a nice feature of the model. As it 
is apparent from Fig.Hl the presence of a dark energy dominated era before the matter 
epoch makes it possible to partially solve the coincidence problem. Indeed, the ratio 
r(z) = Pm(z) I 'pde(z) crosses the unit value two times in the universe history. As a 
consequence, r(z)/r changes by two order of magnitudes (from 0.1 to 10) over the 
redshift range (—0.5, 2.5). Converting the redshift into the dimensionless time r = t/to, 
we see that r(z) stays close to its present day value over most of the universe history 
thus significantly alleviating the coincidence problem. 
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5. Fourth order equivalent theory 

The likelihood analysis performed above has successfully demonstrated that the 
model we have introduced in Sect. 2 is able to reproduce the available astrophysical 
data. This encouraging result makes us confident that our approach to halting the 
eternal acceleration through the introduction of a negative cosmological constant is 
observationally well founded. Nevertheless, a model comprising two dark ingredients 
(the scalar field and the A < term) may be considered as unsatisfactory because 
of the need to theoretically motivate two rather than one new fluids. Moreover, one of 
them is an unusual negative cosmological constant which is difficult to reconcile with 
the classical interpretation in terms of vacuum energy. As such, it is interesting to work 
out a possible reinterpretation of the model in terms of a radically different approach. 

Fourth order theories of gravity have recently attracted a lot of attention as a 
valid alternative to explain cosmic acceleration without the need of any dark energy 
fluid. Moreover, it has been developed a quick method to find out a f(R) theory giving 
the same cosmic dynamics of a given dark energy model. Since all the tests we have 
considered above rely on the Hubble parameter, it is obvious that the reconstructed 
f(R) theory will fit the data in the same way as the model discussed up to now. 

Here, we first summarize the basics of f(R) theories and the method of 
reconstructing the gravity Lagrangian from the Hubble parameter H(z) and then present 
the application to the model we are considering. 

5.1. Basics of f(R) theories 

Much interest has been recently devoted to the so called curvature quintessence 
according to which the universe is filled by pressureless dust matter only and the 
acceleration is the result of the modified Friedmann equations obtained by replacing 
the Ricci scalar curvature R with a generic function f(R) in the gravity Lagrangian. 
The Friedmann equations therefore read [HI [20] : 



3 



Pm 

+ Pc 



(43) 



f'(R) 

2 -)- H lUcurvPcurv > (^^0 

a 

where the prime denotes derivative with respect to R, p curv is the energy density of an 
effective curvature fluid given as : 

1 f 1 

fW) 

and the barotropic factor of the curvature fluid is : 

Rf"(R) + R \Rf"'(R) - Hf"(R) 



Pcurv = TT^y {l lf(R) ~ Rf'(R)] - 3HRf"(R)} , (45) 



w. 



[f(R)-Rf'(R)]/2-3HRf"(R) 



(46) 
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Assuming that there is no interaction between matter and curvature terms, the 
continuity equation for p curv reads [23] : 

, mi \ _ 3H$Q M Rf"(R) _ 3 

Pcurv ~r Oil y 1 ~r W 'curv) f 'curv r „ , , ,3 O \^ ' / 

[J l-^JJ 

which is identically satisfied as can be easily shown using Eq. fH3|) and expressing the 
scalar curvature R as function of the Hubble parameter : 

r=-q(H + 2H 2 ) . (48) 



Combining Eqs.( j43l with Eq.( j44"l) and using the definition of H, one finally gets the 
following master equation for the Hubble parameter [23] : 

H - {3H 2 n M a- 3 + Rf'(R) + R [kf"\R) - Hf\R)] } . (49) 

Inserting Eq. (T4"8"j) into Eq. fT49l . one ends with a fourth order nonlinear differential 
equation for the scale factor a(t) that cannot be analitically solved also for the simplest 
cases (for instance, f(R) oc R n unless dust matter contribution is discarded). Moreover, 
although technically feasible, a numerical solution of Eq. p9l is plagued by the large 
uncertainties on the boundary conditions (i.e., the present day values of the scale factor 
and its derivatives up to the third order) that have to be set to find out a(t) by solving 
Eq.flU]). 

Given these mathematical difficulties, a different approach has been proposed in 
[23] (hereafter CCT) where Eq. p9|) is considered as a way to determine f(R) rather 
than a(t). Rearranging the different terms suitably, CCT obtained a linear third order 
differential equation for / in terms of the redshift z = 1/a — 1 (having set oq = 1) that 
can be easily solved numerically for a given H(z). By this method, it is then possible to 
find out which f(R) theory reproduces the same dynamics of a given dark energy model, 
thus showing a formal equivalence between these two radically different approaches. 

CCT developed the method using the redshift z as integration variable since it is 
common to have an analytical expression for the Hubble parameter as function of z. 
However, this is not the case for the model considered here, so that it turns out to be 
numerically best suited to use t as integration variable. To this aim, we follow [55J to 
rewrite the main formulae in CCT to finally get the equation determining f(t) is therfil: 

^ 3(t) S" + U2 ^% + Hl ^li = - 3H o^MR 2 a' 3 (t) , (50) 

with : 

Hi = 2HR + HR + 2R 2 R~ l - d 3 R/dt 3 , (51) 

U 2 = - (2R + HR) , (52) 

H 3 = R. (53) 

% With an abuse of notation, we write fit) rather than f[R(t)]. 
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where R is given by Eq. (j48j) . In order to integrate Eq. (l501) . we need to specify boundary 
conditions that are more conveniently assigned at the present time. We slightly 
generalize here the discussion presented in CCT. First, let us remember that, in a fourth 
order theory, we may define an effective gravitational constant as G e ff = Gn/ f'(R), with 
Gn the usual Newtonian gravitational constant. Its rate of variation will be given as: 

G eff = 1 f"(R) dR 

G eff t f'(R) dr ■ { ] 

It is quite natural to assume that the effective and the Newtonian gravitational couplings 
take the same values today, so that we get the condition : 

f'(Ro) = 1 . (55) 

Evaluating Eq. (|5^1) for t = t Q (i.e., r = 1), we may determine f"(Ro) provided an 
estimate of (G e // /G e ff) t =t is given. Since in our theory Gn is constant, we may assume 
that the measurements of the variation of Gn [56] actually refers to G e ff and use these 
results to get an estimate of {G e ff/G e ff) t =t - We thus take as our boundary condition: 



having used Eq. (l55|) . Finally, inserting Eqs. fl55l) and (1561) into Eq. (pi5l) and then in (J43] 
evaluated today, we get : 



f(Ro) = QHl ( 1 - n M + f'(R ) 




+ 6H [ — ) f'(R ) . (57) 

t=to 

From Eqs. fl56|) - fl57|) . the following boudary conditions straightforwardly descend: 

^-SL™- (58) 

'fV + ft*) (59) 





t=to \ dt J j._t Q 

that have to be used, together with Eq. (}5T]) . to numerically solve Eq. fl50l) . Combining the 
solution thus obtained for f{t) with R(t) evaluated through Eq. fHSl) . one finally finds 
f{R) thus recovering the higher order theory that mimicks the assigned dark energy 
model. We refer the reader to [221 ESI El] for some interesting examples. 

A preliminary comment is in order here, however. The dark energy model and 
its f(R) duality representation share the same expression for the Hubble parameter 
H(z) (and hence for the scale factor dependence on cosmic time). As such, all 
the observational tests relying directly on H(z) (as the acoustic peak parameter ^4) 
or its integral (as the luminosity distance or the age of the universe) are not able 
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to discriminate between a whatever dark energy model and its f(R) counterpart. 
Nevertheless, such an equivalence only holds at the background level. Indeed, since 
the underlying theory of gravity is radically different, at the perturbative level, the 
two models make clearly distinct predictions. As a consequence, observational tests 
relying on the solution of perturbation equations (such as the growth factor and the 
power spectrum) could be able to break this degeneracy. Actually, at the moment, 
observational determinations of the growth index dkid/dlna are still quite noisy, while 
the uncertainties in the relation betwenn the observationally derived galaxy power 
spectrum and the theoretically predicted matter one still prevent from discriminating 
among dark energy and modified gravity. Moreover, we only rely here on tests depending 
on H(z) so that the f(R) theory we will reconstruct is observationally equivalent to our 
two fluids scenario. 



5.2. Reconstruction of f(R) 

All we need to apply the procedure described above is an expression for the Hubble 
parameter H as function of the cosmic time t. For the model we are considering, this 
may be obtained numerically as described in Sect. 2 so that we only have to choose a 
value for G e ff/G e ff in order to set the initial condition for f"(Ro). We fix G e ff/G e ff = 
in good agreement with most of the estimates in [56]. We have also checked that changing 
this value within the quoted uncertainties does not affect significantly the main results. 

Before presenting the reconstructed gravity Lagrangian, there is a conceptual point 
to clarify. The method described above represents a sort of bridge between two different 
scenarios. In particular, their matter contents could be different, so that we should 
define both Q c ^ rv and fijj^ to denote this quantity in the two different models. In 
principle, there is no reason why fl c ^ rv = VL®p should hold. However, since VL®p is close 
to the fiducial value (VLm — 0.3) suggested by model independent estimates (e.g., from 
galaxy clusters abundance), we take Q c ^/' v = Qffl Should we have chosen a different 
value for Qfj rv , we have had a different reconstructed f(R), but still providing the same 
Hubble parameter as the model with the scalar field and the negative A we have tested 
against data. To be more precise, one could say that our method is not able to recover a 
single f{R), but rather a class of f(R) models parameterized by Q^f v . An astrophysics 
based estimate of Qm is the only way to break this degeneracy making it possible to 
select the most suitable member in this family of fourth order theories. 

With this remark in mind, we have performed the reconstruction of f(R) for the 
model presented in Sect. 2 setting its parameters (Qm,^a,Wo) to their best fit values. 
Rather than reconstructing f(R) directly, it is more instructive to first consider: 

-i 



n 



f(R) 



R 



x 



f(R) 



R 



(60) 

2 = 

which is identically 1 for the Einstein- Hilbert Lagrangian, f(R) = R. Departures of r\ 
from unity therefore quantifies how much the modified Lagrangian should depart from 
the standard one in order to fit the data as well as the input Hubble parameter. 
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Figure 5. rj as function of the redshift z for the f(R) reconstructed from the best fit 
model. Here we define r\ = [f(R)/R] / [f(R)/R] z=0 , see the text for further details. 
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Figure 6. The EoS w curv of the effective curvature fluid for the reconstructed f(R) 
theory. Note that the small fluctuations are only artefacts due to numerical noise. 

To this aim, it is interesting to look at Fig.[5] where we report 77 as function of z 
over the redshift range < z < 6, extending well above the one probed by the data 
we have used. It is apparent that r\ quickly departs from 1, i.e. the reconstructed f{R) 
dramatically differs from the Einsteinian one. This is not unexpected at all. Indeed, in 
order to fit the data with a matter only universe, one should strongly modifies the gravity 
Lagrangian at least for values of R s = R/Ro corresponding to the accelerating expansion 
era. For larger values of z, the universe is decelerating so that one could naively expect 
to recover a standard situation with n approaching 1. This is clearly not the case. 
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Actually, this somewhat counterintuitive result is easily explained reconsidering Fig.H] 
which shows that in our model dark energy dominates again after a finite period of 
matter domination. In order to recover this unusual behaviour, f(R) must still depart 
from the General Relativity Lagrangian. 

Although the reconstructed f(R) is obtained numerically, it is useful to have an 
analytical approximation which we indeed find to be given by : 



with R s = R/Rq, (a, (3,n,m) real parameters to be fitted to the numerical data. We 
have checked that such a formula reproduces the reconstructed f(R) within less than 
4% over the full redshift range we have probed. For the best fit model, we get : 



Changing the values of (Qm, &a, wo) alters these values, but Eq.tlBTT) still provides a very 
good approximation. We stress, however, that the approximating formula flBTl) has been 
tested only over the redshift range (0, 6) so that cannot be extrapolated to larger values. 
As such, one has not to give any weight to the fact that f(R) does not reduce to R in 
the very early universe as expected to recover the BBN succesfull results. 

It is interesting to consider the EoS w curv of the effective curvature fluid which is 
plotted in Fig. [61 This plot shows that the curvature fluid actually behaves as a sort of 
usual cosmological constant with its EoS increasing with resdshift, but still remaining 
quite close to the A value w = —1 over the full redshift range probed by the data. 
This result could be somewhat anticipated noting that the reconstructed f(R) theory 
provide the same H(z) and hence the same AH/H as the starting model. Since this 
quantity stays within 10% over the range < z < 6, it is not surprising that the 
effective curvature fluid behaves as an almost A term. However, we have checked that 
w curv departs more and more from w = — 1 as z increases even if we are not able to 
draw definitive conclusions because of the dramatical increase of the numerical noise. 

As a final comment, it is worth stressing again that the reconstructed f(R) 
provides the same cosmic dynamics (i.e., scale factor and Hubble parameter) of a model 
comprising both a quintessential scalar field and a negative cosmological constant. As a 
consequence, being the background dynamics exactly the same, the f(R) model based 
on Eq. (IBTT) predicts the same behaviour for the luminosity distance as function of the 
redshift z (and hence the same SNela Hubble diagram and gas mass fraction vs z curve) 
and equal numerical values for the acoustic peak and shift parameters. As such, all the 
observational tests in Sect. 3 successfully met by the negative A model are equally 
satisfied by its f(R) duality representation constructed here. Such a fourth order 
theory is therefore able to both fit the available data and avoid eternal acceleration. 
In particular, the net effect is to mimic an unusual negative cosmological constant 
with a modified gravity Lagrangian and no dark energy. This interesting result shows 
that a A < term could result as a net effect of forcing the gravity theory to be an 
Einsteinian rather than a fourth order one. Moreover, halting eternal acceleration is 
therefore possible without adding any ad hoc ingredient to the cosmic pie. 




(61) 



(a,P,n,m) = (-0.38,-0.61,-0.15,0.70) . 
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6. Conclusions 

The unprecedented high quality data accumulated in the recent years have depicted the 
scenario of a spatially flat universe presently undergoing a phase of accelerated expansion 
thus motivating the lot of interest devoted to the search of viable candidates to drive 
this cosmic acceleration. Although many possible models have been put on the ground, 
almost all of them predict a two phase scenario with standard dust matter dominating 
the first epoch of decelerated expansion and a (actual or effective) dark energy fluid 
fueling the present day cosmic speed up. Notwithstanding the underlying mechanism, 
the fate of the universe seems to be yet written. As matter fades away, dark energy 
becomes more and more dominant leading to an eternal accelerationE] so that serious 
problems arise with the formulation of the S matrix and quantum field theory. 

In an attempt to avoid these problems, we have explored here a model made out 
of three rather than two ingredients. Besides the usual dust matter and radiation and 
a scalar field responsible of cosmic acceleration, a negative cosmological constant has 
been added to the energy density budget. While being subdominant during the matter 
dominated decelerated phase and the scalar field dominated accelerated epoch, this new 
component may contrast the action of the quintessential scalar field during the future 
evolution of the universe. Depending on the balance between these two terms, eternal 
acceleration may be avoided leading to a Big Crunch collapse in a finite time. This 
possibility renders the model quite attractive from a theoretical point of view and have 
motivated us to further explore its viability through a comparison with the available 
astrophysical data. To this aim, we have fitted the model to the most recent release of 
SNela data (made out of data coming from the SNLS, ESSENCE and GOODS surveys) 
and the gas mass fraction in galaxy clusters also setting a prior on the acoustic peak 
and shift parameters. This large dataset has been successfully reproduced thus showing 
that the presence of a negative A is fully consistent with the data. As such, our model 
may be considered as a viable solution to the problem of eternal acceleration. 

It is, however, worth investigating whether the same model may be interpreted as 
an effective representation of a radically different scenario. Such a possibility is offered 
by fourth order theories of gravity according to which cosmic acceleration is the result 
of a modified gravity Lagrangian with the scalar curvature R replaced by an analytic 
function f(R) to be determined by the data. Using a method developed by one of us 
[23j EH]) we have reconstructed f(R) by imposing that the cosmic dynamics is the same 
as the one determined by the action of the scalar field and the negative A. Since the 
Hubble parameter is the same in both scenarios, both models are able to fit the data 
in the same way. Moreover, the reconstructed f(R) does not include any cosmological 
constant term thus showing that a negative A could also be the consequence of forcing 
a fourth order Lagrangian to the linear Einsteinian one. 

+ Here, we are not considering phantom dark energy which leads to a superaccelerated expansion 
ending in a Big Rip like singularity. Note, however, that phantom scenarios are affected by the same 
problems of eternal acceleration till the cosmic doomsday. 
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Although encouraging, these results must be considered only preliminary. Indeed, 
the model is affected by some undesirable shortcomings when applied to the early 
universe. First, the model turns out to be matter dominated only for a finite redshift 
range, while its expansion is scalar field driven during most of the high redshift range. 
Moreover, the scalar field EoS converges towards the stiff matter value w = 1 so that 
one could worry about the growth of structure and the CMBR anisotropy spectrum. 
Both these results are not fully unexpected. Indeed, it has yet been demonstrated that 
an exponential potential may lead towards a scalar field dominated universe. Our result 
shows that this conclusion also holds in the case a negative A is present as could be 
predicted in advance noting that this term becomes quickly subdominant at large z 
with respect to the matter and the scalar field ones. One could therefore consider the 
possibility to replace the exponential potential with a different V((f>) using, e.g., the 
sum of two exponential with opposite signs. Suitably weighting the two terms may 
lead to a model which is similar to the present one at small z thus giving the same 
successful match with observations. On the other hand, making the additional term 
dominant for large z may solve the problems of the present model in the early universe. 
We stress, however, that here we have been more interested in looking for a model able 
to halt eternal acceleration. We have therefore paid more attention to the recent and 
future universe so that some problems with the early universe were foreseeable. The 
complicated task of matching the early, late and future expansion looking for a suitable 
potential will be addressed in a forthcoming work. 

Further developments of the approach pursued here are possible from the 
observational and the theoretical point of view. On one hand, we have only tested the 
dynamics of the model, but severe constraints may be obtained investigating structure 
formation (see, e.g., [57] for a discussion of this issue for the expoential potential model). 
To this aim, a first important step may be done solving the perturbation equation in 
the linear regime thus determining the growth factor D = 5/a and the growth index 
/ = din 6/ din a. Since a negative A has the same effect of an attractive gravitational 
potential, it is possible that the collapse of structures is enhanced so that it is worth 
investigating what constraints may be imposed by requiring that this effect does not 
change significantly neither D nor / with respect to the successful ACDM predictions. 
It is worth stressing, at this point, that which are the perturbation equations to solve 
depend on the scenario adopted. Should we interpret the model in terms of f(R) 
theories, the standard theory of linear perturbations must be abandoned and different 
equations have to used. This gives us the possibility to select between the two equivalent 
descriptions retaining the one which predicts D(a) and f(a) most similar to the ACDM 
ones since they are more likely able to fit the large scale structure data. 

On a different ground, one can arguably discriminate between the two components 
dark energy model and its f(R) duality representation based on Solar System tests. 
Indeed, in the low energy limit, the static and spherically symmetric solution of the 
Einstein equations for the two components dark energy model leads to the usual 
Newtonian potential (apart some negligible deviations introduced by the tiny neative A 
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and the scalar field) so that all the classical tests of gravity are automatically verified. On 
the other hand, the low energy limit of f(R) theories may lead to significant deviations 
from the Newtonian potential. As such, these theories are indeed severely constrained 
by the Solar System tests so that some models have been carefully designed just to 
evade this problem [5E]. Although a strong debate is still present on which is the the 
correct way to work out the low energy limit of such theories, it is worth investigating 
this issue with some detail in a future work in order to see whether our approximated 
f{R) expression ( }6Ti) is indeed viable from this point of view. Should this not be the 
case, we could give it away, but still retaining the two components dark energy model 
clS db viable way to halt eternal acceleration. 

A more philosophical comment is in order. Along the work presented here, we have 
first introduced a two components dark energy model comprising a negative A and an 
exponential potential scalar field 0. After successfully testing it against the data, we 
have worked out its fourth order counterpart determining the function f(R) entering 
the modified gravity Lagrangian giving the same cosmic expansion H(z). A guidance 
idea has been the Occam's razor that suggests to give out models with an unnecessary 
large number of unknown elements. The price to pay has been to replace the simple 
general relativistic f(R) = R with a fourth order expression analytically approximated 
as in Eq. flSI]) . One could wonder whether such a complicated Lagrangian may indeed 
be in line with the Occam's razor spirit. Actually, we are at the moment unable to find 
out a fundamental theory leading to our reconstructed f(R), although some Lagrangian 
containing the sum of terms with different powers of R have been studied in literature 
(see, e.g., third reference in [22]). Should a fundamental theory not be found, one can 
indeed consider unnecessary our modification of f(R) and put our theory on the same 
ground as the two components dark energy model we have started from. 

As a conclusive remark, we believe it is worth stressing that halting eternal acceler- 
ation is possible. Using a two component model comprising a quintessential scalar field 
and a negative A represents the easiest way to generate a transient cosmic acceleration. 
Fourth order theories make it possible to interpret both these (somewhat problematic) 
terms as effective manifestations of a different scenario. The reconstructed f(R) may 
then be seen as the missing bridge between the easy world of many components models 
and the mathematically complicated structure of higher order gravity theories. 
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